Pressure and alloying effects on the metal to insulator transition in NiS2- :r Se x . studied 

by infrared spectroscopy 
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The metal to insulator transition in the charge transfer NiSa-xSea; compound has been inves- 
tigated through infrared reflectivity. Measurements performed by applying pressure to pure N1S2 
(lattice contraction) and by Se-alloying (lattice expansion) reveal that in both cases an anoma- 
lous metallic state is obtained. We find that optical results are not compatible with the linear 
Se-alloying vs Pressure scaling relation previously established through transport, thus pointing out 
the substantially different microscopic origin of the two transitions. 

PACS numbers: 71.30.+h, 78.30.-j, 62.50.-p 



Understanding the physics of strongly correlated sys- 
tems is one of the most challenging tasks of condensed 
matter research Besides displaying extremely in- 

teresting physical behavior, their sensitivity to small 
changes in external parameters makes them highly ap- 
pealing for future technological applications. That sen- 
sitivity is attributed to the small value of the electron 
bandwidth in comparison with other relevant energy 
scales as the electron correlation U or the charge trans- 
fer (CT) energy gap. The independent electron approx- 
imation breaks down and materials at half filling can be 
insulators, contrary to the prediction of band theory. 

The cubic pyrite N1S2, which is a CT insulator fol- 
lowing the Zaanen-Sawatsky-Allen classification scheme 
0, is considered, together with vanadium sesquioxide 
V2O3, a text-book example of strongly correlated mate- 
rials. MS2 attracts particular interest as it easily forms a 
solid solution with NiSe2 (NiS2-£ C Se a ;), which, while being 
iso-electronic and iso-structural to N1S2, is nevertheless 
a good metal. A metal to insulator transition (MIT), in- 
duced by Se alloying, is observed at room temperature 
(T) for x « 0.6, and a magnetic phase boundary from 
an antiferromagnetic to a paramagnetic metal is found 
at low T at about x = 1 (see the inset of FigJTJi) [1| . An 
alternative way to induce a metallic state in N1S2 is ap- 
plying a hydrostatic pressure (P). Following Mott's orig- 
inal idea 0] , this technically challenging procedure offers 
the unique opportunity to continuously tune the band- 
width, without introducing impurities or disorder. High- 
P techniques have indeed been used in the past few years 
to investigate the dc transport properties of NiS2- x Sea; 
[3, Q , and a P induced MIT has been observed in pure 
NiS 2 for P > 4 GPa. 



Infrared reflectivity, in particular under pressure, is a 
very suitable probe to address the physics of strongly cor- 
related systems. The investigation of the T-dependent 
optical properties of V2O3 and their theoretical expla- 
nation in terms of coherent and incoherent excitations 
around the Fermi energy (Ej?) represents one of the most 
compelling successes of the dynamical mean field theory 
(DMFT) [a M, p. However, with few remarkable excep- 
tions 0, [l(J Tlllfl^l , infrared investigations of the MIT 
in CT insulators are still rare, and a thorough optical 
study of NiS2- a; Se a ; vs Se-alloying and applied-P is com- 
pletely lacking to the best of our knowledge. In this pa- 
per we fill this gap, presenting room-T reflectivity mea- 
surements over a broad spectral range on 4 compounds 
(x = 0, 0.55, 0.6, 1.2) of the NiS2- £C Se a; series together 
with optical measurements as a function of P on pure 
N1S2. Experimental data are compared with local den- 
sity approximation (LDA) calculations, and the resulting 
scenario for the two MITs is finally depicted. 

The ambient-P nearly normal incidence reflectivity 
R(u) has been measured between 50 and 35000 cm -1 
on well characterized high density pellets of NiS2- a; Se 2 ; 
[l3l [l4| . An in situ evaporation technique was used to 
measure the reference. The high-P study has been per- 
formed using a diamond anvil cell (DAC). A small piece 
of MS2 was loaded inside the gasket hole together with 
KBr as hydrostatic medium. Great care was taken to ob- 
tain a clean sample-diamond interface where reflectivity 
spectra, R s d(u)), have been measured [1 51 ] - The measure- 
ment was performed at the high brightness infrared syn- 
chrotron radiation source SISSI@Elettra (Trieste) jig . 
Further details on the measurement procedures are re- 
ported elsewhere [13, HI] ■ 
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The R(ui) of N1S2 at ambient-P, shown in FigQJ,, is 
nearly flat from 50 to 10000 cm -1 except for weak phonon 
contributions at 260 and 290 cm -1 . On increasing the 
Se-content, R(to) is progressively enhanced at low fre- 
quencies, characteristic of a metallic behavior. The real 
part of the optical conductivity (J\{ui) (FigQ})) has been 
determined through Kramers-Kronig (KK) transforma- 
tions. To this end, standard extrapolation procedures 
were adopted at both high and low frequency [HI, [2(| . 

The optical conductivity of N1S2 is strongly depleted 
at low frequency showing the CT gap (evaluated at Full- 
Width-Half-Maximum of the absorption). This is con- 
sistent with previous optical measurements [3, [2l| at 
about 4000 cm -1 . On increasing the Se-content x, a 
large amount of spectral weight (SW) is transferred from 
high to low frequency through an isosbestic point around 
8000 cm -1 . This suggests the main role played in the 
MIT by electronic correlations [![. As it is better high- 
lighted by the Acri = ai(x) — ai(x = 0) difference spec- 
tra in the inset of FigQJ), the low energy contribution 
to the ai(u>) is made up of two well distinct terms: one 
broad mid-IR band peaked around 2000 cm^ 1 and ex- 
tending up to nearly 8000 cm -1 , and a sharp contribu- 
tion below 500 cm -1 . In analogy with spectra of metal- 
lic V2O3 [y, [7[, the narrow peak is attributed to quasi- 
particle (QP) coherent excitations around Ep, while the 
mid-IR term is associated to optical transitions from the 
QP peak to the upper and lower Hubbard bands. This 
scenario has been confirmed by fitting the <j\ (w) curves 
through a Drude-Lorentz (DL) model. Data can be de- 
scribed by a Drude term plus two Lorentzian oscillators. 
The Drude and the low energy oscillator (centered around 
2000 cm -1 ) describe the coherent and the mid-IR ex- 
citations around Ep, while the remaining oscillator at 
10000 cm -1 mimics the CT and Hubbard transitions. 
The fitting components are reported as thick dashed lines 
in FigQJ) for x = 0.6. 

We turn now to the high-P measurements on MS2. 
The reflectivity at the sample-diamond interface, i? S£ j(w) 
(thick solid lines), is shown in Fig(2k- The strong two- 
phonon diamond absorption provides reliable data only 
above 2000 cm -1 . On increasing the pressure, R s d(u) 
is progressively enhanced at low-frequency showing an 
overdamped behavior (similarly to R^^uj) obtained on 
varying x), as a signature for a correlated bad metal- 
lic state. At high-frequencies all R s d{(^j) converge above 
10000 cm" . In order to evaluate the accuracy of our 
high-P measurements, we compare the R sc i{to) data to 
the expected reflectivity at a sample-diamond interface, 
i?sdVLcalculatedby 

using a procedure previously intro- 
duced OlEl. The calculated R$(w) for NiS 2 (dashed 
lines in FigJ5^) is in good agreement with R s d(co) mea- 
sured in the DAC at the lowest pressure (1.1 GPa) being 
both nearly flat and with a value 20% over the whole 
frequency range. The same procedure has been applied 
to NiS2- 2; Se a ; (x = 0.55, 0.6, 1.2) compounds and the re- 
sulting R°^{oS) are shown in Figl^k for the sake of com- 
parison. We then tried to fit the R s d{u) measured vs P 



1.0 t 



0.8 



„ 0.6 

3, 

2 

0.4 



0.2 
0.0 

3000 

S 

-J 

b 2000 
§ 

1000 



x=0.55 



4-~ 



x=1.2 



=0.6 




1 

x (Se content) 



X=() 



NiS 2 . x Se x 



a) 



3000 

•§• 2000 
b 

T < 1000 

^\x=1.2 
\ x=0.6 


■ 1 1 

1 l^ 1 ^ 


2 4 6xl0 3 
0) (cm" 1 ) _ 


w?--^~x=0.55 ^ 







5000 
0) (cm" 1 ) 



10000 



FIG. 1: (Color online) a) Optical reflectivity of NiS2- a: Sea ; 
for x = 0, 0.55, 0.6, 1.2 at ambient conditions. Inset: 
Phase diagram of NiS2- : rSe :E [lj. Black dots correspond to 
the samples measured in this work, b) Optical conductiv- 
ities from KK transformations. Thick dashed lines repre- 
sent the DL fit of the x = 0.6 sample. Inset: Difference 
Aai = <Ji(x) — <Ji(x = 0) spectra. 



for NIS2 within the same DL framework previously used 
for NiS2- 2; Se 2 ; at ambient-P. Although a certain degree of 
arbitrariness remains in fitting the data over a restricted 
spectral range, a reliable description of the R s d{u) at 
any P is obtained by the sum of a Drude and a mid- 
IR term plus a high-frequency oscillator kept constant at 
all pressures. This fit provides a robust estimate of the 
quasi-particle SW, defined by the sum of the Drude and 
of the mid-IR intensities. This sum remains nearly un- 
changed by varying the fitting parameters over realistic 
ranges. 

The microscopic mechanisms inducing the P- and Se- 
MITs are further investigated by studying the quasi- 
particle SW as a function of the cubic lattice parameter a. 
The lattice is expanded by Se-alloying (22I l23| whereas 
it is compressed by pressure (24|. The x- and the P- 
dependence (up to « 5 GPa) of a have been obtained 
from Ref.[23[ and Ref.[24| respectively. Data at higher 
P have been estimated using the procedure developed in 
Ref . . Through the specific heat results of Ref.[25j, 
which provide a Debye frequency ojd ~ 350 cm -1 well 
comparable with the N1S2 phonon frequencies, we ob- 
tain a sound velocity v s ~ 4300 m/s. As the density 
of MS2 is p = 4455 kg/m 3 , the Bulk modulus results 
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B = p-v 2 s ^ 83 GPa. Assuming B(P) = B + B'P, a{P) 
is finally given by the Birch- Murnagham (B-M) equation 

M 



a{P) = a(0) 



Bn 



*P 



-1/3B' 



(1) 



where B' normally ranges between 4 and 8 [27j . The 
experimental a(P) data, the values estimated from EqQ] 
and those from LDA calculations [28j are in a very good 
agreement as shown in the inset of Fig. [2};. 
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FIG. 2: (Color online) a) R s d(u) measured for NiS2 at high- 
P (thick solid lines) (open symbols show DL fits) and cal- 
culated R^iuj) at sample-diamond interface for NiS2-xSe x 
at selected x (dashed lines). The kink in R s( i(u)) at about 
3700 cm -1 is instrumental. Inset: QP spectral weight vs P 
(bottom) and x (top); top and bottom scales are chosen con- 
sistently with a scaling factor / » 0.14/GPa (see text). Lower 
panels: QP spectral weight (see text) vs lattice constant a for 
N1S2 (b) and NiS2- :c Se :c (c). Green circles are square plasma 
frequency values calculated with LDA (see text). A rescaling 
factor of 1.25 has been used for the comparison with the ex- 
perimental data. The dashed-dotted vertical line marks the a 
value for NiS2 at ambient conditions. Inset: lattice parameter 
vs P. Experimental data from Refl24l (solid circles), calculated 
values using the B-M equation (open circles) and LDA (solid 
triangles). 



The quasi-particle SW, shown for pure MS2 at working 
P in Figl2b and for NiS2- K Se a; at different Se-contents in 
Figl2t, reveals a striking non-monotonic behavior as a 



function of a. Its slow continuous increase for a < 5.57 
A (i.e. at the highest values of P) reflects the progres- 
sive enhancement of the kinetic energy due to the applied 
P and corresponds to a nearly-complete metallization of 
MS2. For a > 5.57 A up to a eq rj 5.68 A (namely the 
lattice parameter corresponding to N1S2 at ambient con- 
ditions), correlation effects get larger and the SW drops 
rapidly to zero as a consequence of the Mott transition. 
On further increasing a above a eq due to the Se-alloying, 
the SW (Fig[2fc) restarts to increase, owing to the onset 
of the Se-induced MIT. 

Despite the opposite behavior of the lattice parame- 
ter vs Se-alloying and pressure, a linear scaling factor 
/ w 0.14/GPa between x and P has been formerly estab- 
lished from low-T dc-resistivity data [J, Q , thus suggest- 
ing an equivalency between the two MITs. However, the 
same / does not apply comparing the optical SW depen- 
dence on P and x. It is indeed clear from the inset of Fig. 
[2^ that the rate of increase of SW is much larger with P 
than with x. The breakdown at finite frequencies of the 
dc linear scaling between x and P suggests that, while a 
metallic state can be obtained from N1S2 both by apply- 
ing P and by alloying Se, this state takes place through 
substantially different microscopic mechanisms, involv- 
ing different redistributions in the electronic density of 
states. A qualitative understanding of the two different 
MITs can be obtained through self-consistent TB-LMTO 
LDA calculations [3(J, see Fig. [3] To this end, we em- 
ployed the iVth order muffin-tin orbital (NMTO) down- 
folding [3l| , and the augmented plane waves plus local or- 
bitals (APW+lo) techniques within the framework of the 
Wien2K code [32|. At ambient P the Ni e g -states with 
a bandwidth W e = 2.1eV and the antibonding ppa*- 
S states are separated by a CT gap (Alda) centered 
around 1.5 eV. Beside this gap, a second LDA-CT gap 
is present between occupied ppn- S states below Ep and 
the e„ 



^-Ni states (see e.g. Ref. [33j). Upon applying pres- 
sure, the lattice contracts and the entire band-structure 
around Ep is renormalized; e.g. by a factor of 1.13 at 
P = 10 GPa (dashed curve in Fig. [3h,). All features of 
the DOS stay the same, in particular, the energy scales, 
(W e and Alda, are rescaled by a factor 1.13. Hence, 
the bandwidth-gap ratio W e / Alda remains nearly con- 
stant. On the other hand, the interaction U can be as- 
sumed to be constant, so that W e JU increases by the 
factor 1.13, triggering a bandwidth-controlled MIT (BC- 
MIT). In the case of Se-substitution, the lattice expands 
(instead of shrinking) due to the larger atomic radius of 
Se ions. This leads to a very complementary scenario: 
While the changes of the e ff -bandwidth W e are negligi- 
ble, the CT gaps shrink (see Fig. [5}}). Assuming U again 
to be constant, the driving force for the MIT is now the 
the reduction of the charge transfer gap (W Bg / Alda in- 
creases, as W eg remains basically unaffected). 

Let us now turn back to the optical experiment in Fig. 
[2 Deep inside the metallic phase, we can expect cor- 
relations to be weak and LDA to give the proper an- 
swer. Within LDA we calculated the "square plasma 



4 



frequency", which is defined as the average over the FS 
of the squared velocities. In Fig. [2b-c we compare the 
square plasma frequency to the experimental SW. Be- 
cause of the different changes in the bandstructure under 
pressure and upon Se alloying, LDA gives non-monotonic 
behavior and very different slopes in agreement with ex- 
periment, see Fig. HJd-c). For the insulating N1S2 com- 
pound and close to the phase transition, electronic corre- 
lations are not negligible and we performed LDA+DMFT 
[34| calculations. To this end, the NMTO bandstructure 
was downfolded to effective Ni e 9 -states and the corre- 
lations in these two orbitals were treated by means of 
DMFT. For U > 3 J (U being the intra-orbital Coulomb 
interaction between the Ni e 9 -states and J the Hund cou- 
pling) , the two e 9 -orbitals split and a gap opens for N1S2 . 
This insulating LDA+DMFT solution results in a very 
strong suppression of the square plasma frequency (as 
well as of the SW) as indicated by the vertical arrow in 
Fig. db-c. 
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FIG. 3: (Color online): a) Total LDA-DOS for NiS 2 . Solid 
line represent the DOS at ambient P, dashed lines at 10 GPa. 
b) Total LDA-DOS for the NiSe2 compound. In the case of 
NiS2, the 10 GPa DOS can be rescaled by a factor 0.88 on 
that at ambient-P, resulting in an increase of W £g /U, keeping 
M4 b /Alda fixed. In contrast, Se-substitution results in a 
decrease of W eg /U and an increase of W Bg / Alda due to the 
shrinking of the charge-transfer gap Alda. 



Besides important similarities between the P- and Se- 
dependent phase diagrams [1, the present optical 
study reveals that a simple linear scaling between P and 
x , as that indicated by transport, does not hold at finite 
frequencies. This suggests that the two MITs rely on dis- 
tint microscopic mechanisms. These mechanisms can be 
understood theoretically in terms of the two fundamental 
parameters for the MIT in a CT insulator: Under pres- 
sure, W eg /ALDA — const, and W eg /U increases, trigger- 
ing the MIT; in contrast upon alloying Se, the increase of 
W eg / Alda is responsible for the MIT, whereas W e /U 
even decreases. This makes NiS2- a: Se 3 ; under pressure an 
ideal system for the study of the MIT in a CT strongly 
correlated system. 
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